Root-mean-square deviation (bioinformatics)

The root-mean-square deviation (RMSD) is the measure of the average distance between the atoms (usually the backbone atoms) of superimposed proteins. In the study of globular protein conformations, one customarily measures the similarity in three-dimensional structure by the RMSD of the Cα atomic coordinates after optimal rigid body superposition.

When a dynamical system fluctuates about some well-defined average position the RMSD from the average over time can be referred to as the RMSF or root-mean-square fluctuation. The size of this fluctuation can be measured, for example using Mössbauer spectroscopy or nuclear magnetic resonance, and can provide important physical information. The Lindemann index is a method of placing the RMSF in the context of the parameters of the system.

A widely used way to compare the structures of biomolecules or solid bodies is to translate and rotate one structure with respect to the other to minimize the RMSD. Coutsias, et al. presented a simple derivation, based on quaternions, for the optimal solid body transformation (rotation-translation) that minimizes the RMSD between two sets of vectors.[1] They proved that the quaternion method is equivalent to the well-known Kabsch algorithm.[2]

Contents

The equation

\mathrm{RMSD}=\sqrt{\frac{1}{N}\sum_{i=1}^N\delta_i^2}

where δ is the distance between N pairs of equivalent atoms (usually Cα and sometimes C,N,O,Cβ).

Normally a rigid superposition which minimizes the RMSD is performed, and this minimum is returned. Given two sets of n points \mathbf{v} and \mathbf{w}, the RMSD is defined as follows:

\begin{align} \mathrm{RMSD}(\mathbf{v}, \mathbf{w}) & = \sqrt{\frac{1}{n}\sum_{i=1}^{n} \|v_i - w_i\|^2} \\ & = \sqrt{\frac{1}{n}\sum_{i=1}^{n} ({v_i}_x - {w_i}_x)^2 %2B ({v_i}_y - {w_i}_y)^2 %2B ({v_i}_z - {w_i}_z)^2} \end{align}

An RMSD value is expressed in length units. The most commonly used unit in structural biology is the Ångström (Å) which is equal to 10–10m.

Uses

Typically RMSD is used to make a quantitative comparison between the structure of a partially folded protein and the structure of the native state. For example, the CASP protein structure prediction competition uses RMSD as one of its assessments of how well a submitted structure matches the native state.

Also some scientists who study protein folding simulations use RMSD as a reaction coordinate to quantify where the protein is between the folded state and the unfolded state.

See also

References

  1. ^ Coutsias EA, Seok C, Dill KA (2004). "Using quaternions to calculate RMSD". J Comput Chem 25 (15): 1849–1857. doi:10.1002/jcc.20110. PMID 15376254. 
  2. ^ a b Kabsch W (1976). "A solution for the best rotation to relate two sets of vectors". Acta Crystallographica 32 (5): 922–923. doi:10.1107/S0567739476001873. 

Further reading

External links